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ABSTRACT 

We present the redshift distribution of a complete, unbiased sample of 24 fim sources down to 
/„(24/xm) = 300 /^Jy (5a). The sample consists of 591 sources detected in the Bootes field of the 
NO AO Deep Wide-Field Survey. We have obtained optical spectroscopic redshifts for 421 sources 
(71%). These have a redshift distribution peaking at z ~ 0.3, with a possible additional peak at 
z ~ 0.9, and objects detected out to z — 4.5. The spectra of the remaining 170 (29%) exhibit no strong 
emission lines from which to determine a redshift. We develop an algorithm to estimate the redshift 
distribution of these sources, based on the assumption that they have emission lines but that these 
lines are not observable due to the limited wavelength coverage of our spectroscopic observations. The 
redshift distribution derived from all 591 sources exhibits an additional peak of extremely luminous 
(Lg_iooo//m > 3 x 10 12 L Q ) objects at z ~ 2, consisting primarily of sources without observable 
emission lines. We use optical line diagnostics and IRAC colors to estimate that 55% of the sources 
within this peak are AGN-dominated. We compare our results to published models of the evolution 
of infrared luminous galaxies. The models which best reproduce our observations predict a large 
population of star- formation dominated ULIRGs at z > 1.5 rather than the AGN-dominated sources 
we observe. 

Subject headings: galaxies: distances and redshifts — galaxies: evolution — galaxies: formation — 
infrared: galaxies 



1. INTRODUCTION 

An important goal in extragalactic astronomy is to un- 
derstand when, where, and how the stellar content of the 
universe formed. Based on rest-frame ultraviolet and op- 
tical measures, the star formation rate per comoving vol- 
ume is approximately constant from z ~ 6 to z rs 2 
and then smoothl y decreases by a factor of ten until 
the present epoch dGallego et alj|1995t iLilly et al.lll996t 
Madau et all Il99ft IConnollv et all 119971: rTrever et all 
1998 flFlores et al.lfl 999: Stcidel et al. 19991: IWilson et all 
20021 : iGiavalisco et alj 12004). There is also increasing 
evidence that galaxies with large stellar masses formed 
most of their stars earlier than lower-mass galaxies 
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(iTinslevI 1 19681 ICowie et all 119961 : iHeavens et all l2004t 
Uuneau et al.ll2005D . 

Observations of the far-infrared background radia- 
tion with the Cosmic Background Explorer ( COBE ; 
Pueet et al lll996HFixsen et alj|1998t lHauser et al.lll998f 
Dwek et alj|1998t lLagache et al.lll999HLagache fc Puge 
200tt iFinkbeiner et alJl2000f ) show that the infrared en 



ergy density is comparable to the combined UV, visible, 
and near-infrared energy density, indicating that half of 
the light originally emitted in the UV and optical has 
been interce pted by dust, reproc essed, and re-emitted in 
the infrared (|Gispert et al.| [2000). What kinds of galax- 
ies contribute to the infrared peak in the extragalactic 
background light? The earliest mid-infrared observa- 
tions of extragalactic sources revealed a class of galaxies 
which emit the bul k of th e ir luminosity in th e infrar ed 
(Low & Klcinmann 1968; Kleinmann fc Low! fL970al rhl: 
iRieke fc Low! fl972) . Early observations with the In- 
frared Astronomical Satellite (IRAS) showed that lumi- 
nous infrared galaxies (LIRGs; 10 11 < L 8 _i O oo^m(L0) < 
10 12 ) and ultraluminous infrared galaxies (ULIRGs; 
Ls-iooo/im > 10 12 L ) do not contribute substantially 
to the total infrared emission in t he local (z < 0.05) uni- 
verse (jSoifer fc Neugebauenfl99l[h 

However, there is strong evidence that infrared-bright 
populations become more significant at higher redshift. 
This e vidence includes nu mber counts obtained with 
IRAS (lHacking et alj 119871) and the Infrared Space Ob- 
servatory - (ISO; lOliver et alJ liap; lAltieri etafl 119991: 
l Aussel et all 119991: lElbaz et all 119991 : iGruppioni et all 
l2002t : Metcalfe et alj [20031 ). Spectroscopic follow-up of 
these sources has directly confirmed the increasing sig- 
nificance of the infrared-bright population out to z ~ 0.3 



2 



Desai et al. 



(|Kim et all 119981: ISerieant et all l2001h . a redshift set 
by the sensitivity limits of the existing mid-infrared 
surveys. Progress in constraining the population of 
infrared-luminous galaxies at higher redshi ft has been 
made via surveys at submillimeter (e . g. ISmail et aTl 
119971 IChapmanet al.l [2001L iBlain et all 120021 ) and ra- 
dio (ICowie et al . 2004) wavelengths. These data pro- 
vide further evidence of strong evolution in the infrared- 
luminous population, but the number of sources probed 
are small, and may suffer from selection biases. For in- 
stance, spectroscopy of submillimeter sources has relied 
on accurate positions obtained from radio counterparts, 
biasing spectroscopic samples to source s that are bright 
in bo th the submillimeter and the radio (|Chapman et all 
2005). The sensitive Multiband Imagi ng Photometer for 
the Spitzer Space Telescope (MIPS; iRieke et all 12004 
I Werner et al.l [20041 ) has allowed the detection of many 
faint infrared sources. Number counts at 24 /im deter- 
mined by Spitzer observations confirm the strong evolu- 
tion previously observed in the infrared luminous popula- 
tion (ICharv et all 120041: iDole et all 12004 iMarleau et~aTl 
120041: iPapovich et al.ll200l ~ 

Although number counts provide valuable informa- 
tion on the evolution of the infrared population, they 
do not allow a robust distinction between luminosity 
and density evolution. For this, redshifts are needed. 
The redshift distribution is shaped by the evolution of 
galaxy SEDs and luminosity functions, which in turn de- 
pend on the physics of galaxy evolution: star formation, 
AGN activity, dust obscuration, mergers, and feedback 
(|Lagache et al.ll2~005|) . The redshift distribution is there- 
fore a useful constraint on the importance of these pro- 
cesses. 

Several attempts have been made to estimate the red- 
shift distribution of 24 /im sources using photometric 
redshifts. For sources with J„(24/im) > 83 /iJy, these 
studies consistently reveal a primary peak in the red- 
shift distribution at z < 1, with small differences in 
the precise location of this peak attributable to cosmic 
variance ([P erez- Gonza lez etall 120051: iLe Floc'h et~aT1 
120051: ICaputi et all 12001 IPapovich et all l2007t ). Such 
a peak can be attributed to strong luminosity 
and/or density e volution in the infra r ed population 
out to z ~ 1 (ICharv fc Elbazl 120011: lLagache et all 
12004 iGruppioni et al.ll2005t iPerez-Gonzalez et al.ll2005f ). 
However, contradictory results are found for the ex- 
istence of a predicted secondary peak at z ~ 2 
(jPerez-Gonzalez et al.1 120051 : ICaputi et all 120061 ). Such 
a secondary peak may be expected from a significant 
high-redshift population of either highly absorbed AGN- 
dominated sources or PAH-rich starburst-dominated 
sources. These would have a local peak in their SEDs 
at «8 /im, due either to continuum adjacent to deep ab- 
sorption features or to line emission, respectively. This 8 
/im peak would enhance the detectability of these sources 
at z « 2 in 24 /im surveys. Uncertainties in the pho- 
tometric redshifts may be responsible for the contra- 
dictory observational results regarding the existence of 
this peak, since previous studies have few (<25) spec- 
troscopic redshifts above z = 1.5. Photometric esti- 
mates of the redshift distribution for more shallow sur- 
veys (/„(24/im) > 200 / iJy) also produce mixed re- 
sults for the z ~ 2 p eak (|Rowan- Robinson et al.l [20051 
iBabbedge et al.l [20061 ). Recent spectroscopic surveys of 



24 /im sources are complete only to about the same red- 
shift (z ~ 1) reliably probed by ph otometric redshifts 
(|Papovich et al.1 120061 : fMarleau et al.ll2007t) . 

In this work, we present new measurements of the red- 
shift distribution of 24 /im sources down to 300 /iJy. 
Our sample is selected from 23 noncontiguous areas 
spread across 8.74 deg 2 of the 24 /im map of the Bootes 
field of the NO AO D eep Wide-Field Survey (NDWFS; 
iJannuzi fc Devi Il999f ) . Our results are therefore mini- 
mally affected by variations in the density of the 24 /im 
source population due to structures on the scale of ~0.5 
deg 2 . Although our sample is modest (591 sources; com- 
parable in size to the smaller studies mentioned above), 
a larger fraction (71%) have spectroscopic redshifts. The 
remaining 29% exhibit no spectral features from which to 
determine a redshift. Given the limitations of photomet- 
ric redshifts for 24 /im sources at z > 1, we have devel- 
oped an algorithm for estimating the redshift distribution 
of this population from their continuum spectra. Briefly, 
this described (algorithm in detail in 35. 2[) is based on 
the assumptions that bright 24 /im sources are under- 
going strong star formation and/or AGN activity, and 
the emission lines resulting from this activity are able to 
penetrate the surrounding dust. These assumptions hold 
true for «70% of the 24 /im sources in our sample, and 
also hold fo r very obscured local ultraluminous infrared 
sources (e.g. [Sanders et aT1ll988atlKleinmann et afll988 |: 
Armus et al.l Il989t lAllen et all 119911: ICutri et all Il994 
Veilleux et al.lll995l T7997t IDuc et al.lll997l ). Under these 
assumptions, the lack of observed emission lines for cer- 
tain sources is due to the limited wavelength coverage of 
our spectroscopic observations. 

We describe our targeting strategy in our targeting 
completeness in Sj21 and our spectroscopic observations 
in 52) Our results are presented in SJSJ In particular, we 
present the redshift distribution of emission-line galaxies 
( 35. ip . the redshift distribution of galaxies that exhibit 
no emission lines ( 35. 2|) , and the redshift distribution for 
the entire sample ( 35. 3p . In fj6j we discuss the shape of 
the overall redshift distribution and compare our obser- 
vations to existing models of the evolution of infrared 
galaxies. Finally, we summarize our conclusions in fj7[ 
All calculations are carried out assuming a cosmology 
with H = 70 km s" 1 Mpc" 1 , Q, Q = 0.3, and ft A = 0.7. ' 

2. TARGETING STRATEGY 

We have conducted an optical spectroscopic survey of 
~600 24 /im sources in order to determine their redshift 
distribution. The spectroscopic targets were drawn from 
a catalog of ~22,000 sources detected in a MIPS survey 
of 8.74 deg 2 of the Bootes field of the NDWFS. The 5a 
point-source depth of the imaging is 300 /iJy, and this 
sets the flux limit of both the detection catalog and our 
spectroscopic survey. 

The spectroscopy was carried out using the Deep Imag- 
ing M ulti-Object Spectrograph (DEIMOS; iFaber et all 
2003) and the Low R esolution Imaging Spectrometer 
(LRIS; lOke et al.l[i~995l ). both on the W. M. Keck IO- 
meter telescopes. In total, we observed 11 DEIMOS and 
12 LRIS masks. Mask positions were chosen to maximize 
the number of high-priority targets, which are defined to 
be sources with /„(24/im) > 750 /iJy and faint optical 
magnitudes (R > 24 mag or f v {R) < 0.74 /iJy). Figure[Tj 
shows the positions of the masks on the sky. The remain- 
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ing targets on each mask were selected from the sources 
with /^(24/im) > 300 jitJy, without regard to R - [24] 
color. To /„(24/zm) > 300 /iJy, 24 /zm sources have a sur- 
face density of about 2400 deg -2 . There are ~50 per 0.02 
deg 2 DEIMOS mask and «25 per 0.01 deg 2 LRIS mask, 
resulting in a total of 818 over the 0.33 deg 2 combined 
mask area. We targeted 544, or 67%. (Our final sample 
includes 47 supplementary sources described in Section 
[31) The untargeted sources were excluded mainly due to 
slit collisions. In the case of slit collisions, we preferen- 
tially targeted faint (but not invisible) optical sources. 
Although their combined area is modest, our spectro- 
scopic masks span 8.74 deg 2 of the Bootes field. This 
makes our survey less subject to the effects of cosmic 
variance due to structures on the scale of approximately 
half a square degree. 

In designing the spectroscopic masks, we attempted to 
position the slits on the optical (i?-band) counterparts of 
the 24 /zm targets. For reference, the NDWFS i?-band 
imaging has a 3c point-source depth of ^26 Vega mag- 
nitudes. 13 The position of the optical counterpart was 
determined in the following way. For each 24 /im source, 
the SExtractor i?-band catalog was searched for nearby 
sources. For extended sources, we match detections in 
different bands when their centroids are within an ellipse 
defined by the second-order moments of the optical light 
distribution (as defined by SExtractor). For both com- 
pact and extended sources we match detections in dif- 
ferent bands when their centroids are within 1" of each 
other. When matching bright sources, where centroids in 
different bands can be offset or optical sources are well 
resolved, matching MIPS sources to SExtractor catalogs 
has some advantages over simple i?-band aperture pho- 
tometry at the position of the MIPS source. When their 
are multiple i?-band detections which meet our match- 
ing criteria, our code opts for the closest centroid. The 
plots in this paper show magnitudes computed within a 
4" diameter aperture. In some cases, we find no i?-band 
counterpart. After visual inspection of these cases, we 
perform i?-band aperture photometry at the MIPS po- 
sition. The optical positions determined by the above 
method were used as a starting point for designing the 
slit masks. The final slit positions were adjusted based 
on visual inspection of the i?-b and and 24 p,m images , 
as well as additional imaging (Eisenha rdt et alj 120041 ) 
acquired with Sp itzer's Infrared Array Camera (IRAC; 
iFazio et alJl200l . The IRAC images were useful in iden- 
tifying ambiguous optical counterparts, since the MIPS 
targets are all detected by IRAC, which has a much 
smaller PSF than MIPS. In addition, sometimes faint 
optical counterparts can be identified by eye even if they 
did not appear in the i?-band SExtractor catalog. 

3. TARGETING COMPLETENESS 

In the following, the term "Keck spectroscopic tar- 
gets" refers to both the high-priority targets and the flux- 
limited targets described above. The inclusion of high- 
priority sources as well as the preference for faint i?-band 
sources in cases of slit collisions mean that our sample 
of Keck spectroscopic targets is not completely unbiased. 
Figure [2] shows the differential and cumulative distribu- 

13 See http:/ /www. noao.edu/noao/noaodeep/ for more informa- 
tion regarding the depth and coverage ot the NDWFS. 



tions of R- [24] colors for all of the /^(24/zm) > 300 /zJy 
sources within the Keck survey region and for the subset 
of these that were targeted for spectroscopy. A two-sided 
Kolmogorov-Smirnov (KS) test shows that the 24 /im 
flux density distributions of these two sets of galaxies 
are indistinguishable. In contrast, their i?-band magni- 
tude and R — [24] color distributions are different at a 
statistically significant level, in the sense that our target 
list includes an excess of optically faint, red sources. 

To correct this bias, we supplement our Keck spectro- 
scopic targets with 47 optically bright (R < 23 mag) 
sources within the Keck survey region that were not tar- 
geted with LRIS or DEIMOS, but which have redshifts 
from the AGN and Galaxy Evolution Survey (AGES; 
Kochanek et al. in preparation). AGES is a spec- 
troscopic survey of galaxies (I < 20 mag) and AGN 
(I < 21.5 mag) in the NDWF S Bootes field carried out 
with the Hectospec instrument (jFabricant et al.|[2005[ ) on 
the MMT. We selected the supplementary AGES sam- 
ple randomly under the requirement that each 0.1 R- 
band magnitude bin at R < 23 mag have the same sam- 
pling fraction as the average Keck sampling fraction at 
R > 23 mag (74%). The differential and cumulative 
distributions of R — [24] colors for the combined Keck 
and supplementary AGES targets are shown in Figure [2j 
Two-sided Kolmogorov-Smirnov tests confirm that the 
Keck+AGES sample is fairly drawn from the population 
of sources with /„(24/im) > 300 /iJy that lie within the 
Keck spectroscopic survey area, in terms of / 1/ (24/xm), R- 
band magnitude, and R — [24] color. Thus we conclude 
that the Keck+AGES targets are an unbiased sample of 
sources selected at 24 /*m. 

4. SPECTROSCOPIC OBSERVATIONS 

The DEIMOS and LRIS optical spectroscopic followup 
of 24 /jm sources was carried out in May 2004 and May 
2005. The DEIMOS observations for each mask are com- 
posed of 3-4 frames lasting 20 or 30 minutes each, result- 
ing in total exposure times of 90-120 minutes per mask. 
The 600 line mm -1 (7500 A blaze) grating was used with 
the GG400 blocking filter, resulting in a mean spectral 
dispersion of 0.65 A pixel -1 . The central wavelength 
was set to 7500 A. We typically used 1.0" wide by 10" 
long slitlets, corresponding to a FWHM spectral resolu- 
tion of 4.6 A. Although the wavelength coverage varies 
with slit position, a typical DEIMOS spectrum ranges 
from 5200-10200 A. We used the DEEP2 data reduction 
pipeline 14 to perform cosmic ray removal, flat-fielding, 
co-addition, sky-subtraction, and wavelength calibration. 
The one-dimensional s pectra were extracted using stan- 
dard IRAF 15 routines (|Todvlll993l ) . 

The LRIS observations are also composed of 3-4 frames 
lasting 20 or 30 minutes each, yielding total exposure 
times of 60-120 minutes per mask. One mask was ob- 
served for a total of only 20 minutes. The telescope was 
offset by a few arcseconds between frames. The slitlets 
were 1" wide by 11" long. LRIS was configured with 
the 5600 A dichroic. On the blue side, we used the 400 
line mm -1 grism with a central wavelength of 3400 A, 

14 http:/ /alamoana. keck. hawaii.edu/inst/deimos/pipeline. html 

15 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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resulting in a spectral dispersion of 1.09 A pixel -1 and 
a FWHM spectral resolution of 6.9 A. On the red side, 
we used the 400 line mm -1 grating with a central wave- 
length of 8500 A, resulting in a spectral dispersion of 
1.86 A pixel -1 and a FWHM spectral resolution of 7.5 
A. Our LRIS spectra typically have a wavelength cov- 
erage of 2200-8200 A. We used an IRAF package called 
BOGUS 16 to reduce the two-dimensional spectra. The 
one-dimensional spectra were extracted using standard 
IRAF routines. 

Relative spectrophotometric calibration for both the 
DEIMOS and LRIS observations was perform e d us- 
ing observations of Wolf 1346 fMassev etall 119881 : 
iMassev fc Gronwalllll990l:lOkdll990l) . ~ 

5. RESULTS 

Our primary observational result is the redshift dis- 
tribution of 24 /mi sources down to 300 /iJy. In the 
following three subsections, we present the redshift dis- 
tributions of the sources with detected emission lines, of 
the sources without detected emission lines, and of the 
entire sample. 

5.1. Redshift Distribution of Galaxies with Detected 
Emission Lines 

Each of the one- and two-dimensional DEIMOS and 
LRIS spectra were visually inspected for spectral lines 
from which redshifts could be measured. Some example 
spectra are shown in Figure [3j and many more will be 
provided in a forthcoming paper (Soifer et al. in prepa- 
ration). Out of the 544 sources that were targeted, 368 
spectra (69%) yielded redshifts. Two or more spectral 
lines were typically available to determine a redshift. 
In 44 cases (12%), only one spectral line was observed. 
Eight of these had absorption features in the Keck spec- 
trum or additional features from AGES spectra that 
helped to determine the redshift. In the remaining spec- 
tra, the single observed emission line was identified based 
on its shape and width as Lya in eight cases, MgIIA2798 
in three cases, and [OIIJA3727 in 25 cases. These single- 
line sources are included in our sample with spectroscopic 
redshifts. 

Of the 176 galaxies that we targeted with either 
DEIMOS or LRIS but for which we were unable to de- 
termine a redshift, six of these were bright enough in the 
optical (R = 21.16, 16.58, 19.42, 22.15, 20.12, 18.30 mag) 
to have redshifts from the AGES Survey. The AGES 
redshifts are z = 1.799, 0.076, 2.230, 1.509, 0.545, and 
1.069. The DEIMOS spectrum of the z = 0.076 source 
did not yield a redshift because it was saturated. The 
LRIS spectrum of the z = 0.545 source showed absorp- 
tion lines, but identifiable emission lines were lost in the 
gap between the red and blue sides. The differing wave- 
length coverage of the Keck spectrographs compared to 
Hectospec explains the remaining sources. The AGES 
data bring the number of Keck targets with redshifts 
up to 374. Hereafter, we refer to the redshifts of all 
Keck targets with either Keck or AGES spectroscopic 
redshifts as the Keck sample. In contrast, the Keck tar- 
gets supplemented with sources with bright optical fluxes 
from AGES, as described in will be referred to as the 
Keck+AGES sample. 

16 https:/ /zwolfkinder.jpl. nasa.gov/~stern/homepage/bogus. html 



Figure [4] shows the i?-band magnitudes versus 24 /im 
flux densities for the Keck+AGES targets, for the sources 
in the spectroscopic survey region that were not included 
in the Keck+AGES target list, for the Keck+AGES tar- 
gets with spectroscopic redshifts, and for the Keck tar- 
gets without spectroscopic redshifts. This figure illus- 
trates that the 24 /im flux density is only weakly cor- 
related with the i?-band magnitude, and that the Keck 
targets for which we were unable to determine spectro- 
scopic redshifts are mostly faint in the optical (R > 21 
Vega mag), but span a range in 24 /jm flux densities. 

Figure[S]shows the success rate in determining redshifts 
for Keck and Keck+AGES targets as functions of 24 /im 
flux density, i?-band magnitude, and R— [24] color. This 
success rate is a strong function of i?-band magnitude, 
and drops below 50% for R > 23.5 mag. This dependence 
on i?-band magnitude propogates into a dependence on 
the R — [24] color. In contrast, the redshift success rate 
is a much weaker function of 24 /im flux density. 

Figure [S] shows the 24 /im flux density, i?-band mag- 
nitude, and R — [24] color versus redshift for the sources 
with spectroscopic redshifts. The 24 /im flux density is 
not tightly correlated with redshift. In contrast, both the 
i?-band magnitudes and the R — [24] colors of the major- 
ity of these same sources correlate strongly with redshift 
out to z = 1. At larger redshifts, the scatter increases 
substantially. 

For each source with an emission line redshift, Fig- 
ure [7] shows the vL v luminosity at the rest wavelength 
probed by the MIPS 24 /tm filter at the observed red- 
shift. The luminosity- dependent SEDs developed for the 
iCharv fc Elbazl (|200l models are used to indicate the 
regions of this plot occupied by normal galaxies, LIRGs, 
ULIRGs, and HyperLIRGs. This figure shows that the 
24 /im sources in our sample span a large range in bolo- 
metric luminosities. 

5.2. Redshift Distribution of Galaxies without Emission 

Lines 

We were unable to determine emission line redshifts for 
170 Keck targets (124 DEIMOS targets and 46 LRIS tar- 
gets). Visual inspection of the optical, IRAC, and MIPS 
imaging suggests that the optical counterpart may have 
been incorrectly targeted in up to 20% of these sources. 
Generally, this uncertainty is due to the faintness of the 
optical counterpart combined with the large 24 /im beam 
size. A quarter of the time, the uncertainty stems from 
a choice between multiple optical counterparts. Despite 
this uncertainty, we use the following algorithm to esti- 
mate the redshift distribution of these sources. 

First, we assume that these sources lie at z < 4.5. This 
is the highest emission-line redshift found in our survey. 
This redshift cut-off is further justified by the fact that 
a z > 4.5 source with /„(24/mr) > 300 /iJy would have a 
luminosity o f Lg-iooo/ini > 3 .5 x 10 13 L Q (assuming the 
templates of lCharv fc E lbaz (2001)). Galaxies with such 
extreme luminosities should be exceedingly rare. 

Second, we assume that 24 /im sources that are well- 
detected in the bluest NDWFS filter (B w < 25 mag) 
must lie at z < 3. This ensures that some part of the Bw 
filter (which covers the approximate wavelength range 
3500-4750 A) is sampling rest-frame wavelengths long- 
ward of the Lyman break. 
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Third, we assume that all bright MIPS sources should 
be undergoing strong star formation and/or AGN ac- 
tivity, and should therefore harbor the conditions re- 
quired to produce strong emission lines, including Lya, 
[OIIJA3727, H/3, and Ha. Although red sources that 
are bright at 24 /im must be dusty, these emission 
lines are able to penetrate the dust for the major- 
ity of our targets. This is likely due to the fact 
that the dust obscuration is patchy, allowing opti- 
cal/UV emission lines to be observed in very ob- 
scured sources. In fact, many very obscured local ul- 
traluminous infrared galaxies, including Arp 220, dis- 
play strong emission lines which dominate their opti- 
cal spectra (IDuc et all 119971; IVeilleux et al.lll997l 1 19951: 
Cutri et all 11994 lAllen et all 119911: lArmus et al.lll989l ; 
Kleinmann et alj|1988t iSanders et all ll988bT T We as- 
sume that this is the case for the targets without detected 
emission lines as well, and that the lack of observed emis- 
sion lines is due to insufficient wavelength coverage rather 
than to intrinsically weak or completely obscured emis- 
sion lines. Below we describe how we rule out those red- 
shifts for which any one of Lya, [OII]A3727, H/3, or Ha 
could have been detected. 

We assume that a line could have been detected at a 
particular wavelength (corresponding to a specific red- 
shift) if the limiting rest-frame equivalent width at that 
wavelength is smaller than a fiducial rest-frame equiv- 
alent width, which we take to be 10 A for each of the 
four emission lines we are considering. For each spec- 
trum, we determined the limiting rest-frame equivalent 
wi dth as a function of wavelength following Equation 4 
in lHogg et all (|1998l ): 

1 /2 

EQW lim , rost (A) = f^j [i + z(A)]" 1/2 . (1) 

Here, r\ is the minimum signal-to-noise ratio required for 
the line to be detected, and was taken to be three for 
both LRIS and DEIMOS spectra, and for all four emis- 
sion lines. The variable r(A) is the signal-to-noise ratio 
of the extracted spectrum as a function of wavelength. 
It was computed by convolving the extracted spectrum 
with a Gaussian kernel characterized by a FWHM of 200 
A, and dividing the result by the noise within a mov- 
ing bin of 13 pixels in the (unconvolved) extracted spec- 
trum. The spectral dispersion, Ai, is 0.65 A pixel -1 for 
DEIMOS, 1.09 A pixel -1 for the blue side of LRIS, and 
1.86 A pixel" 1 for the red side of LRIS. The parame- 
ter AA is the rest-frame FWHM of the emission line in 
question, which we take to be 10 A for all four emission 
lines. Finally, z is the redshift which corresponds to the 
wavelength at which we are computing the rest-frame 
limiting EQW, and is given by z = (A/A rcs t) — 1, where 
A rcst = 1216,3727,4861,6563 A for Lya, [Oil], H/3, and 
Ha, respectively. 

After ruling out those redshifts at which these emission 
lines could have been detected for each individual source, 
we constructed a redshift probability distribution func- 
tion for the ensemble by summing the number of sources 
that could lie at a given redshift, and normalizing the 
result to a total probability of unity. This is based on 
the assumption that each galaxy has a uniform proba- 
bility of lying at any redshift that is not ruled out by 
the above three arguments (that all 24 /xm sources lie at 



z < 4, that all 24 yum sources that are well-detected in 
Bw must lie at z < 3, and that all 24 /im sources must 
have strong emission lines that have gone undetected due 
to limited spectral coverage) . Figure [8] shows the result 
for all of the targets without redshifts, and for those ob- 
served with DEIMOS and LRIS separately. Based on the 
lack of the [OII]A3727 emission line in their spectra, most 
of the sources observed with DEIMOS are unlikely to lie 
at z < 1.5. In addition, few lie at z > 3.5, or else we 
would have observed Lya. In contrast, the redshift dis- 
tribution of the sources observed with LRIS is closer to 
uniform. The similarity of the top and bottom panels of 
Figure[5]reflects the fact that most of the targets without 
spectroscopic redshifts were observed with DEIMOS. 

Figure [8] does not take into account the fact that the 
mean luminosity of galaxies in a flux-limited survey in- 
creases with redshift and the luminosity function de- 
creases rapidly at the high-luminosity end. Qualitatively, 
this means that a galaxy does not have an equal prob- 
ability of lying at every redshift not ruled out by the 
above algorithm. A given galaxy is more likely to be an 
L* galaxy at moderate redshift than five times more lu- 
minous at a higher redshift. However, given the limited 
photometric information available to derive a luminos- 
ity at a common wavelength for each galaxy, we do not 
attempt to quantify this qualitative statement. 

5.3. Redshift Distribution 

The redshift distribution of the Keck+AGES targets 
is shown in Figure [9l The contributions of the following 
categories of sources are indicated separately: the Keck 
targets with emission-line redshifts; the supplementary 
AGES sample selected to correct the bias of the Keck 
spectroscopic targeting against bright optical sources; 
and Keck targets without spectroscopic redshifts, assum- 
ing that they all have emission lines that could not be 
observed due to the limited wavelength coverage of the 
Keck spectrographs and that they have an equal proba- 
bility of lying at any redshift in the range < z < 4.5 
that is not excluded by the lack of detected emission lines 
(see Section 15.21) . Figure [9] shows that the galaxies with 
spectroscopic redshifts peak at z ~ 0.3. The majority of 
the targets without emission lines cannot lie at this red- 
shift or [OII]A3727 would have been observed. Instead, 
most of these sources are constrained to lie at 2 < z < 3, 
where [OII]A3727 lies longward of the observed spectrum 
and Lya lies shortward of it. These sources appear in the 
redshift distribution as a weak second peak at z ~ 2. 

6. DISCUSSION 
6.1. Shape of the Redshift Distribution 

The redshift distribution of 24 /zm sources is expected 
to be shaped both by evolution in the infrared luminosity 
function and by major features in the mid-infrared spec- 
tra of luminous infrared galaxies. These in turn depend 
upon fundamental processes governing galaxy evolution, 
such as star formation, AGN activity, dust obscuration, 
galaxy mergers, and feedback from supernovae and AGN. 
Thus, the redshift distribution is an important constraint 
on models of these processes. 

In the absence of any evolution, the number of sources 
observed in a given redshift bin depends only on the vol- 
ume probed by that bin, the range of luminosities probed 
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by that bin, and the shape of the z = luminosity func- 
tion. The first two effects act in opposite directions: vol- 
ume increases with redshift, leading to an increase in 
the observed number of galaxies with redshift. However, 
the flux limit means that a decreasing range of luminosi- 
ties, and therefore a decreasing number of galaxies, is 
observed with increasing redshift. How does an evolv- 
ing luminosity function affect the redshift distribution? 
A decrease in either the space density or the character- 
istic luminosity of infrared sources with redshift would 
cause fewer galaxies to be observed at high redshift. Con- 
versely, an increase in the space density or characteristic 
luminosity would appear as a rise in the redshift distri- 
bution. 

In addition to a strong dependence on the evolving lu- 
minosity function, the observed redshift distribution is 
also affected by prominent features in the mid-infrared 
spectra of infrared sources. A strong emission feature 
can increase the observed broadband MIPS 24 /xm flux 
density of a galaxy. Thus, objects with strong emission 
features can be pushed into a flux-limited sample, even 
if they would not be targeted based on their continuua 
alone. Conversely, sources with strong absorption fea- 
tures can fall out of flux-limited samples. The strongest 
mid-infrared spectral features in LIRGs and ULIRGs are 
the silicate absorption feature at rest-frame 9.7 /im and 
the Polycyclic Aromatic Hydrocarbon (PAH) emission 
features at rest-frame wavelengths 3.3, 6.2, 7.7, 8.6, 11.3, 
12.7, 16.3, and 17 /xm. The redshift distribution of 24 /jm 
sources may display troughs and peaks at redshifts cor- 
responding to the presence of these absorption and emis- 
sion features in the observed 24 /xm bandpass, and there- 
fore contains information about the types of galaxies at 
each redshift. In this subsection, we consider the most 
prominent features of the observed redshift distribution 
shown in Figure [9l compare to the results of previous 
observations, and discuss their possible origins. 

Figure [5] shows a strong increase in the number of 
sources fr om z = to z ~ 0.3 . A ri se is also seen at 
z < 1 bv iPerez-Gonzalez et aD (|2005f ). iLe Floc'h et alJ 
(2005), and ICaputi et al] ([2006). However, the peak in 
these studies occurs at higher redshifts (z ~ 0.8), pre- 
sumably because they have a deeper 24 /im flux density 
limit (83 /xJy versus our 300 /xJy). As discussed above, 
a rise in the redshift distribution could be caused by a 
strong spectral feature entering the 24 /xm bandpass, by 
evolution in the luminosity function, or simply by the 
fact that larger volumes are probed at higher redshifts. 
There are no strong spectral features that could be caus- 
ing the z ~ 0.3 peak. The detailed modeling that is 
required to distinguish between the remaining scenarios 
is beyond the scope of this paper. However, in ^6.21 we 
compare our observations to existing phcnomcnological 
models of the evolving luminosity function. 

There is potentially another peak in the redshift distri- 
bution at z ~ 0.9. However, it is very weak, and has not 
been seen in the previous measurements of the redshift 
distribution described in [JTJ Given the small areas of 
the previous surve ys, this could very well be due to cos- 
mic variance (see ISomerville et al.1 (|2004[ )). Additional 
possibilities include the fact that most previous efforts 
were either incomplete at these redshifts, or based on 
photometric redshifts, which would tend to smooth out 
features in the intrinsic redshift distribution. If real, this 



peak could be the result of the 12.7 /im PAH feature and 
the 12.8 /im [Nell] emission line passing through the 24 
/xm bandpass. 

Between z ^ 0.8 and z ~ 1.6, the redshift distribution 
ex hibits a steep decline, whic h ha s also been reporte d 
bv lPerez-Gonzalez et al.l (|2005| ) and lCaputi et al.l((2006h . 
Such a decline could be the result of a flattening or a drop 
in either the density or the characteristic luminosity of 
infrared galaxies. A decline could also result if a signif- 
icant fraction of sources with / i ,(24/xm) > 300 /iJy have 
deep silicate absorption features at rest-frame 9.7 /xm. 
Some highly absorbed sources may fall below the sam- 
ple's flux limit as the absorption feature enters the 24 
/im bandpass. This would result in a broad dip centered 
at z ~ 1.5, the redshift at which the center of the broad 
absorption feature coincides with the center of the wide 
bandpass. 

Also visible in the redshift distribution is a possible 
peak at z ~ 2. However, Figure [9] clearly shows that this 
peak is composed mostly of sources for which we were un- 
able to determine spectroscopic redshifts due to the lack 
of observed emission lines. We remind the reader that 
we constrained the redshift distribution of such sources 
by assuming that 1) these non-detections were due to the 
limited wavelength coverage of optical spectrographs; 2) 
24 /im sources that are well-detected in the bluest ND- 
WFS filter (Bw < 25 mag) must lie at z < 3 so that 
some part of the B\y filter is sampling rest-frame wave- 
lengths longward of the Lyman break; and 3) there are 
no 24 /xm sources at redshifts larger than z — 4.5 (see 
q5.2[) . The exact location of the z ~ 2 peak is uncertain 
because we have assumed that a given galaxy is equally 
likely to lie at any redshift not ruled out by the previous 
three assumptions. 

In the future, it will be possible to provide stronger 
constraints on the redshift distribution of the mid- 
infrared sources which lack spectroscopic redshifts by in- 
corporating priors associated with the mid-infrared lumi- 
nosities, optical luminosities, and colors of the sources. 
However, adding priors to our estimate at the present 
time is complicated for many reasons. First, there is no 
well-measured luminosity function for 24 /xm sources at 
high redshift. The current best luminosity functions are 
estimated using pho tometric redshifts a nd do not extend 
beyond z ~ 2 (e.g. ICaputi et ail 120071 hereafter C07). 
Second, the sources for which we are attempting to esti- 
mate the redshift distribution are the reddest ones whose 
SEDs are poorly fit by galaxy templates and for which 
photometric redshifts are the least reliable. Third, the 
optical luminosity distribution of 24 /xm sources is not 
known as a function of redshift. Although it is expected 
that fainter sources will lie at higher redshift (Figure 0, 
it is not clear how to implement this expectation and 
what the contribution to the counts might be, if any, 
from low-redshift dwarfs. Finally, we note that our spec- 
troscopic redshift sample shows a clear trend of R — [24] 
color with redshift, but the relation is broad and it is not 
clear if the sources that lack spectroscopic redshift mea- 
surements partake in the trend or are outliers. We made 
an attempt to implement a simple luminosity prior by ex- 
trapolating the evolving luminosity function as reported 
by C07. The results for all of the Keck+AGES targets as 
well as for the subset of these with spectroscopic redshifts 
are shown in Figure [9] Using the ILe Floc'h et al.l £2005) 
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luminosity function extrapolated beyond 2 = 1 produces 
similar results. Compared to the C07 luminosity func- 
tion, it more closely reproduces our spectroscopic redshift 
distribution at z < 1 , with the exception of the small dip 
at z ~ 0.75. As one might expect, a simple luminosity 
prior results in weighting the redshift probability distri- 
bution of a given source to lower redshift. This would 
be balanced by the optical brightness prior, which would 
push faint sources to higher redshift. Depending on the 
redshift evolution assumed, the strength of the bump at 
z ~ 2 can be reduced in strength. However, since there 
is no knowledge of the forms or weights of the priors dis- 
cussed, we present the redshift distribution in Figure [9] 
as the measured distribution (white histogram) and our 
simple prediction for the sources without spectroscopic 
redshifts as estimated from the spectroscopic data with 
no priors. These should be interpreted as likely bracket- 
ing the true redshift distribution of the 24 /im sources. 

The presence of a peak at z ~ 2 has been difficult to 
establish in previous studies of th e 24 /im redshift dis- 
tribut ion. It was not detected by iPerez^G onzalc z et al.l 
(l2005h but i t was detected with high significance by 
iCaputi et al.l (|2006f ). Uncertainties in the photometric 
redshifts may be resp o nsible for th is difference. N either 
iPerez-Gonzalez et all (|2005l) nor ICaputi et all (|2006l ) 
have very many (<25) spectroscopic redshifts above z = 
1.5 (See Figures 13 and 14 in the former and Figures 1 
and 5 in the latter). As a result, the scatter in the pho- 
tometric redshifts at z > 1.5 has not been quantified. 
Without this information, an assessment of the reality 
of a peak at z ~ 2 based on photometric redshifts is 
difficult. (For comparison, our Keck spectroscopic sur- 
vey includes spectroscopic redshifts for 46 galaxies at 
z > 1.5.) Both P05 and C05 focused on sources down to 
83 fj,3y. Our study suggests that even sources brighter 
than 300 /iJy display this peak, although this is based on 
assumptions regarding the nature of the sources without 
detected emission lines. 

Neural net and template-fitting techni ques have been 
used to calculate photometric redshifts (jBrodwin et alJ 
l2006h for NDWFS sources detected by IRAC, including 
the sources presented in this work. These photometric 
redshifts are highly accurate for IRAC-selected galaxies 
and AGN out to z — 1.5 and 2 = 3, respectively. An 
analysis of 24 /im selected galaxies with spectroscopic 
redshifts shows that their photometric redshifts are accu- 
rate out to z ~ 1, where the 24 /im sample is is dominated 
by normal galaxies and LIRGs. At higher redshifts, cor- 
responding to ULIRGs, the photometric redshifts be- 
come unreliable (the redshift dispersion in a /(I + z) is 
0.13 with 95% clipping). Since it is likely that the Keck 
sources without spectroscopic redshifts lie at z > 1, we 
cannot use these photometric redshifts to determine their 
redshift distribution. 

There are two reasonable interpretations of the z ~ 2 
peak in the redshift distribution. The first is that there 
is a significant population of infrared-luminous sources 
with strong 7.7 and 8.6 /im PAH features. These sources 
have increasing detectability at 1.6 < z < 2 and de- 
creasing detectability at 2 < z < 3 as the PAH emis- 
sion shifts in and out of the 24 /im bandpass. The 
sources comprising the z ~ 2 peak have /„(24/im) > 300 
/iJy, and are therefore ULIRGs with infrared luminosi- 
ties Lg-iooo/iin > 3 x 10 12 L Q (assuming the tem- 



plates of ICharv fc Elbazl (|2001ft ). In the local universe, 
galaxies with these luminosities typic a lly have low PAH 
equivalent width s (|Lutz et all 119981 : iTran et all 120011 : 
iDesai et"al"1l2007l ). However, there are sources at z ~ 2 
that are known to be both extremely luminous and have 
strong PAH features dLutz et all 120051: lYan et all 120051 : 
iMenendez-Delmestre et alll2007t ISaiina et al [ 12001 7 in- 
dicating that the relationship between PAH equivalent 
width and luminosity may be evolving in the sense that 
extremely luminous galaxies at high redshift have spec- 
tral prop erties similar to le s s luminous galaxie s at low 
redshift (|Sajina et all l2007t IDesai et al J 120071 ). Such 
an evolution might be expected if the PAH equivalent 
width is a function of various evolving sources properties 
(e.g. the distribution of star formation and the presence 
of an AGN). If the relationship between luminosity and 
spectral properties does evolve, then it is reasonable that 
a significant fraction of the sources composing the z ~ 2 
peak have large PAH equivalent widths. This is an im- 
portant effect to include in phenomenological models of 
galaxy evolution, such as those discussed in the next sub- 
section. 

A second explanation for the shape of the redshift dis- 
tribution at z > 1.6 is that there is a large population 
of sources at these redshifts with deep 9.7 /im silicate 
absorption features. As discussed above, similar sources 
could be at least partly responsible for the decline in 
the redshift distribution between z ~ 0.8 and z ~ 1.6 
because of decreasing detectability as the 24 /im band- 
pass samples the absorption feature. In contrast, these 
sources are more likely to be detected at z ~ 2, where 
the 24 /im bandpass is sampling rest-frame wavelengths 
shortward of the absorption feature. In this scenario, 
the z ~ 2 peak represents the reappearance of highly ab- 
sorbed sources as the silicate absorption feature moves 
out of the 24 /im bandpass. The decrease at z > 2 would 
then be due to a combination of a decline in the flux 
density at wavelengths shorter than 8 /im and the fact 
that fewer galaxies meet the flux density cut at higher 
redshift in the absence of evolution. 

We have qualitatively argued that the z ~ 2 bump 
could consist of some combination of sources with large 
PAH equivalent widths and sources with power-law con- 
tinuua accompanied by significant silicate absorption. Is 
there any evidence that either of these populations domi- 
nates? Figure [9] shows that a small fraction of the sources 
in the z ~ 2 peak have been spectroscopically confirmed. 
Of the 33 Keck targets with 1.6 < z < 3, 73% show spec- 
troscopic signatures of AGN activity (lines broader than 
2000 km s _1 in the rest frame or the presence of high- 
ionization lines such as NVA1240, CIVA1549, HeIIA1640, 
[NeV]A3346, or [NeV]A3426). The remaining 27% are 
of unknown type. Of the AGN, 50% are Type 1 AGN 
(with lines broader than 2000 km s _1 ). Thus, a signifi- 
cant fraction of the spectroscopically confirmed sources 
in the z ~ 2 peak are AGN. However, these represent a 
small fraction of all of the galaxies in the peak, leaving 
the nature of the majority unknown. 

IRAC colors have been used by several authors to iden- 
tify sources that have power-law S EDs and therefore 
presumably contain po werful AGN (jStern et all 120051 : 
lLacv et all I2004L l2007t ). Using the photometry from 
the IRAC Shallow Surv ey of the NDWFS Bootes field 
(jEisenhardt et all l2004f ). Figure [TO] shows the IRAC 
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color-color diagram for all of the Ke c k targ ets without 
spectroscopic redshifts. iSajina et al.1 (|2005h have mod- 
elled the IRAC color-color diagram, and have delineated 
regions expected to be occupied by starbursts and AGN 
at various redshifts. Given the uncertainties involved, 
the redshift distribution implied by the IRAC colors is 
comparable to that deduced from the spectral coverage 
of the spectra. For example, the fraction of starforming 
galaxies at z < 1.6 suggested by the IRAC colors is sim- 
ilar to the total number of "no-z" sources expected to lie 
at z < 1.6 according to Figure [8j 

Given this concurrence, we assume that the vast ma- 
jority of sources in the IRAC AGN wedge lie at z > 1.6. 
According to Figure [5J 75% of sources at z > 1.6 are 
expected to lie in the z ~ 2 peak, which we define to 
include all galaxies in the range 1.6 < z < 3. If the 
fraction of AGN-dominated galaxies at z > 1.6 reflects 
the fraction within the z ~ 2 peak, then about 63% x 
75% = 47% of sources without spectroscopic redshifts 
that lie within the z ~ 2 peak (50 out of 106 sources) are 
AGN-dominated. 

We have already mentioned that 73% of the 33 spectro- 
scopically confirmed sources within the z ~ 2 peak show 
optical spectroscopic signatures of AGN activity. Figure 
[TU] shows that 85% of these fall within the IRAC AGN 
wedge. In fact, 82% (27 of 33) of all of the spectroscopi- 
cally confirmed sources within the z ~ 2 peak lie within 
the AGN wedge. Only two AGES sources lie within the 
z ~ 2 peak (see Figure |9]). 

Based on the proportion of sources within the AGN 
wedge, we estimate that (50 + 27) / (106 + 33) = 55% of 
sources (including those with and without detected emis- 
sion lines) within the z ~ 2 peak are AGN-dominated. 
As mentioned in Section 1, Caputi et al. (2006) have 
used photometric redshifts to estimate the redshift dis- 
tribution of 24 /im sources down to 80 fiJy. As in the 
current work, they also detect a peak at z ~ 2, which 
they interpret as the result of a population of starburst- 
dominated sources. This is not necessarily in contradic- 
tion with our results, since their flux limit is significantly 
fainter than ours. As shown in Brand et al. 2006 and 
Dey et al. 2008, the fraction of sources dominated by 
star formation is expected to increase with decreasing 24 
jum flux density. Indeed, Figure 5 in Caputi et al. 2006 
shows that most of the sources in their z ~ 2 peak fall 
below the 24 fim flux density limit o f our survey. 

The Infrared Spectrograph (IRS; IHouck et all l200l 
on board the Spitzer Space Telescope has been used 
to obtain mid-infrared spectra of a sizeable number 
of optically faint (R > 23 ma g) , infrared bright 
(,f i y(24^tm) > 750 /iJy) sources (IHouck et al.l 2005: 



Yan et aLpOOllWeedman et alJl2006tlSaiina etall2007 : 
Dev et aill2008l ). The IRS spectra of such sources tend 



to have power-law continuua indicative of AGN activity, 
and some of these also exhibit silicate absorption fea- 
tures implying z ~ 2. Near-infrared spectra of a similar 
sample also reveal high redshifts and rest-f rame optical 
spect roscopic signatures of AGN activity ([Brand et al.1 
2007). At least some of the sources in the z ~ 2 peak 
have similar i?-band magnitudes and 24 /iin flux densi- 
ties as these samples (see shaded histograms in Figure 
[5]). However, many of the sources in the z ~ 2 peak are 
substantially fainter. There is some evidence that the 
contribution of star formation to the bolomctric lumi- 



nosity of a galaxy increase s with decreasing 24 fim flux 
density ([Brand et all 120061 : iDev et all [20081 ). If this is 
the case, then it is possible that the fainter sources in the 
z ~ 2 peak display strong PAH emission, which is typical 
in starburst galaxies. Near-infrared or mid- infrared IRS 
spectroscopy of the Keck targets without redshifts would 
help clarify the presence of the z ~ 2 peak and the en- 
ergy generation mechanisms of the sources that comprise 
it. 

6.2. Comparison with Models 

In the previous subsection, we qualitatively discussed 
the possible origins of the observed features of the red- 
shift distribution shown in Figure [HI We found that 
quantitative models are required for a robust interpre- 
tation of these features. In this subsection, we compare 
our observations to the predictions of several such models 
which have been found in the literature. 

In Figure [TTJ we compare the observed redshift distri- 
bution to that predicted by four models of the evolution 
of mid-in frared sources. The models are descr i bed i n 
detail by iCharv fc Elbazl (l200ll). iLagadie et all ([200l . 
iGruppioni et al.1 (|2005l ). and iPearsonl ( 20051 ): hereafter 
referred to as CE01, L04, G05, and P05, respectively. 
Below, we briefly describe the major elements of these 
models. 

CE01 evolve the local mid-infrared luminosity function 
to fit the spectrum of the cosmic infrared background and 
the galaxy counts at mid-infrared, far-infrared, and sub- 
millimeter wavelengths. They constructed a family of 
luminosity-dependent template SEDs spanning the lu- 
minosities of normal galaxies, starbursts, LIRGs, and 
ULIRGs. In their model, all galaxies with Ls-iooo/zm > 
10 102 L Q evolve, while only 5% of galaxies with lower 
luminosities evolve. The luminosity evolution goes as 
(1 + z) 4 5 up to z = 0.8, after which it remains constant. 
The density evolution goes as (1 + z) 1 ' 5 out to z = 0.8, 
after which it goes as (1 + z)~ 0A . 

L04 follow the evolution of "normal" and starburst 
galaxies. Each of these galaxy types is represented by 
luminosity-dependent SEDs. The magnitude of the lu- 
minosity and density evolution in these populations was 
chosen to fit the existing observed mid- and far-infrared 
number counts, the mid- and far- infrared redshift dis- 
tributions, the far-infrared luminosity functions, the far- 
infrared background, and its fluctuations. The resulting 
model features a very high rate of evolution of the lumi- 
nosity density of starburst galaxies, peaking at z ~ 0.7 
and remaining constant up to z = 4. The normal galax- 
ies evolve out to z = 0.4, after which their luminosity 
density remains constant. 

G05 start with the 15 /im luminosity function of galax- 
ies and AGN. They assume that four populations con- 
tribute: "normal" galaxies, starburst galaxies, Type 1 
AGN, and Type 2 AGN. They constrain the luminosity 
and density evolution of these populations by requiring 
a match to the observed 15 ^m number counts. In this 
model, the normal population does not evolve. Mean- 
while, the starburst population evolves strongly out to 
z = 1, as (1 + z) 3,5 in luminosity and as (1 + z) 3,8 in 
density. Type I AGN evolve in luminosity as (1 + z) 2 - 6 
out to z = 2, with no evolution at higher redshifts. Type 
II AGN evolve in luminosity as rapidly as (1 + z) 2,6 out 
to z = 2 and do not evolve at higher redshifts. 
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P05 present two phenomenological models tracking the 
evolution of four classes of galaxies: normal, starburst, 
LIRG/ULIRG, and AGN. Both models are anchored to 
observed, type-dependent, local mid-infrared luminosity 
functions. In the "bright-end" model, the starbursts and 
AGN evolve in both luminosity and density as (1 + z) 3 - 3 
out to z = 1, and do not evolve at higher redshifts. 
The LIRG/ULIRG component evolves in luminosity as 
(1 + z) 2 - 5 and in density as (1 + z) 3 5 out to z = 1, 
and does not evolve at higher redshifts. In the "burst" 
model, the starburst population behaves as it did in the 
"bright-end" model. However, the LIRG/ULIRG com- 
ponent evolves in density as (1+z) 7 out to z = 1 and then 
more slowly in both density and luminosity to higher red- 
shifts. 

For a fair comparison with the observations, all of the 
models shown in Figure [TT1 have been computed down to 
/^(24/im) = 300 /iJy and have been normalized to the 
number of 24 /im sources targeted by our Keck observa- 
tions (591). The observed redshift distributions shown 
in Figure [11] have been binned in redshift identically to 
the models for easy comparison. Because the models all 
have different binsizes, which are also different from the 
binsize used in Figure O the observations look slightly 
different in each panel and in Figure [9] 

In comparing the models to the observations in Figure 
rm we note that the detailed structure of the observed 
redshift distribution at z > 1.5 is not well-constrained, 
for the reasons described in §5.21 However, the relative 
number of sources at z < 1 and z > 1 should be a robust 
statistic with which to compare models and observations, 
since the vast majority of the Keck sources without spec- 
troscopic redshifts are extremely unlikely to lie at z < 1. 
The observations indicate that 45% of the Keck+AGES 
targets lie at z > 1. In the CE01 and G05 models, 
only ^25% of the sources lie at z > 1. The templates 
used in these models may be insufficient to reproduce 
the observations. The CE01 model lacks both highly 
absorbed AGN-dominated sources and high-luminosity 
sources with strong PAH features, both of which could 
induce a peak at z ~ 2, as discussed in tj6.ll The G05 
model includes no dependence of PAH strength on lumi- 
nosity. Although they include AGN-dominated sources, 
none have deep silicate absorption. 

In contrast, the L04 and P05 models correctly predict 
the fraction of sources at z > 1. In both models, these 
high-redshift sources are starburst-dominated ULIRGs 
with strong PAH features. However, we suggested in §6.11 
that absorbed power-law (AGN-dominated) SEDs likely 
also contribute to a trough at z ~ 1.5 which appears 
as a peak at z ~ 2. Although the relative numbers of 
z < 1 and z > 1 galaxies in the L04 and P05 models 
are consistent with our observations, further modelling is 
needed to determine whether an equally consistent model 
could be built with a higher fraction of AGN-dominated 
sources at z > 1. In addition, the exact position of the 
z ~ 2 peak is critical to testing the L04 and P05 models. 
This would require additional spectroscopy in the near- 
or mid-infrared. 

7. CONCLUSIONS 

We have conducted a redshift survey of 591 bright 
(/„(24/zm) > 300 /iJy) MIPS sources in the NDWFS 
Bootes field. We obtained emission line redshifts for 71% 



of our targets, and have developed an algorithm to derive 
the redshift distribution of the remaining 29%. Using the 
entire distribution, we conclude: 

1. The redshift distribution of bright 24 /im sources 
peaks at z ~ 0.3. There are no strong spectral fea- 
tures entering the 24 /im bandpass at the appro- 
priate redshifts to cause this peak. It is therefore 
likely caused either by evolution in the luminos- 
ity function or by the fact that larger volumes are 
probed at higher redshifts, possibilities which can 
be tested with models. 

2. There is a marginally-detected additional peak at 
z ~ 0.9. If real, it could be due to the fact that 
24 /im surveys are more sensitive at this redshift to 
LIRGs with strong 12.7 /im PAH features and/or 
12.8 /im [Nell] line emission. Additional data are 
needed to test the validity of this peak. The clean- 
est method would be to continue our spectroscopic 
survey of 24 /im sources to build up the sample 
size. Another approach would be to compute pho- 
tometric redshifts for all of the 24 /im sources in 
the Bootes field, trading redshift accuracy (reason- 
able for photometric redshifts less than z ~ 1) for 
a larger sample size. 

3. We find weak evidence for another peak at z ~ 2, 
consisting mainly of sources for which we were 
unable to determine emission- line redshifts, but 
whose redshift distribution we have estimated by 
assuming that the lack of emission lines is due to 
the limited wavelength coverage of the Keck spec- 
trographs. Such a peak would be expected due to 
the increased sensitivity to starforming ULIRGs at 
redshifts where their 7.7 and 8.6 /im PAH features 
pass through the 24 /zm bandpass. Alternatively, 
the peak could be the result of decreased sensitivity 
to heavily obscured AGN at redshifts where their 
9.7 /zm silicate absorption features pass through the 
24 /im filter. Based on the arguments presented in 
the next paragraph, we suggest that the latter is 
an important effect. 

4. We suggest that 55% of the sources in the z ~ 2 
peak are AGN-dominated, and the remaining 45% 
are starburst-dominated. This suggestion is based 
mainly on the diversity of IRAC colors among the 
Keck targets without redshifts. In addition, AGN 
signatures are present in the optical spectra of 73% 
of the few (33) spectroscopically-confirmed sources 
within the z ~ 2 peak. Finally, sources similar 
to the brightest (/^(24/im) > 750 /zJy) Keck tar- 
gets without spectroscopic redshifts have AGN-like 
near- and mid-infrared spectra. 

5. Of the four existing phenomenological models of 
galaxy formation that we considered, those that 
include a significant number of PAH-rich (and 
therefore presumably star-formation dominated) 
ULIRGs at z > 1.5 best reproduce the relative 
number of galaxies observed at z < 1 and z > 1. 
However, it remains to be shown that models in- 
cluding a significant fraction of AGN-dominated 
ULIRGs at z > 1 cannot also fit the observations. 
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Fig. 1. — Schematic of Keck spectroscopic survey. Points indicate the locations of all Bootes sources with /„(24^tm) > 300 /^Jy, the 
sample from which our optical spectroscopic targets were drawn. Overlaid are the positions of the LRIS and DEIMOS masks that comprise 
the spectroscopic survey. Mask positions were chosen to target rare sources with the most extreme infrared-to-optical flux density ratios. 
The final redshift distribution is corrected for this selection as described in J2] The total area of the spectroscopic survey is 0.33 deg 2 . 
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Fig. 2. — Differential (top) and cumulative (bottom) R — [24] color distributions. Each of the distributions is normalized by the number of 
sources within the represented sample. Because of the targeting selection described in J2] the Keck targets include a disproportionately small 
fraction of sources with blue R. — [24] colors. Kolmogorov-Smirnov (KS) tests show that supplementing the Keck targets with additional 
redshifts from AGES results in a sample that is fairly drawn from the full 24 fim Bootes population in terms of R — [24] color. KS tests 
also confirm that the 24 fim flux density and R-band magnitude distributions of the Keck+AGES sample are consistent with those of the 
full Bootes population. 
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The fourth row shows a DEIMOS spectrum with only one emission line, but multiple absorption features allow for a robust redshift 
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Fig. 7. — Luminosity (uL v ) versus rest-frame wavelength for Keck+AGES targets with spectroscopic redshifts. Plotted quantities 
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Fig. 8. — Estimated redshift probability density function (black line) and Poisson errors (grey shaded region) for Keck targets without 
spectroscopic redshif ts. T hese functions were calculated by ruling out redshifts based on the non-detections of Lya, [OIIJA3727, H/3, and 
Ha, as described in £|5.2I This method predicts that the redshift distribution of the "no-z" sources peaks between 1.5 < z < 3. 
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Fig. 9. — Observed redshift distribution of 24 fim sources down to /„(24) = 300 /ijy. The redshift distribution peaks at z ~ 0.3, 
with a marginally-detected additional peak at z ~ 0.9. We find weak evi denc e for another peak at z ~ 2, based on estimates of the 
redshift distribution of the Keck targets without spectroscopic redshifts (see i|5.2l l. The redshift distribution implied by the C07 Lg_ iooo/mi 
luminosity function is also shown, normalized to the number of sources in the entire Keck+AGES target list, and also to the subset with 
spectroscopic redshifts. The C07 luminosity function implies that the sources without spectroscopic redshifts lie at z < 1. We argue in £|5.3I 
that this is unlikely. 
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Fig. 11. — Observed redshift distribution of 24 /im sources compared to the CE01, L04, G05, and P05 models. The models and their 
associated Poisson errors are shown as shaded histograms. P05 present two variants of their model: "burst" and "bright end" . All models 
have been normalized to the total number of Keck+AGES targets. The highest-redshift bin shows the number of counts at rcdshifts greater 
than its leftmost extent. The observed redshift distribution, binned identically to the models, is shown as points with lrr Poisson error bars. 
None of the error bars include any contribution from cosmic variance. The L04 and P05 models, both of which predict a large fraction of 
PAH-rich ULIRGs at z > 1.5, best reproduce the observations. 



